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OBJECT 

To  investigate  the  use  of  high-speed  computer  techniques  for  the 
analysis  of  aircraft  personnel  catapult  performance. 


SUMMARY 

The  high-speed  computer  has  been  used  in  the  analysis  and  predic- 
tion of  catapult  performance.  By  mean.*  of  an  iterative  numerical  inte- 
gration procedure,  simultaneous  records  of  thrust-time,  ve locity-t ime , 
and  displacement -time  can  be  generated  for  a given  catapult  system  in 
a matter  of  a few  minutes.  Preliminary  studies  indicate  that  catapult 
performance  with  combinations  of  high  frictional  drag  (such  as  rail 
friction)  and  high  acceleration  fields  (caused  by  gravity  and  aircraft 
maneuvers)  can  be  evaluated  with  these  new  techniques. 


AUTHOR  I ZAT  I ON 


00  471.4/2750,  FA  471.4/5971,  24  Aug  45 


INTRODUCTION 


This  report  derives  a technique  for  using  a high-speed  digital  computer  to  obtain 
the  ballistic  solution  for  an  aircraft  personnel  catapult.  Thrust,  velocity,  and 
stroke  are  given  as  functions  of  time.  This  technique  may  also  be  used  for  similar 
cartridge  actuated  devices.  The  procedure  consists  of  using  the  equations  of  constraint 
of  the  system,  i.  e.,  the  equations  of  motion,  the  energy  equation,  the  equation  of 
state  or  gas  law,  and  the  equation  of  burning,  together  with  the  initial  conditions, 
in  a stepwise,  iterative,  numerical  integration  procedure.  In  the  solution  thus  gener- 
ated, the  following  assumptions  are  made: 


1.  Gas  leakage  and  initial  friction  are  negligible. 

2.  A linear  burning  equation  is  used. 

3.  The  equation  of  state  for  an  ideal  gas  is  used. 

4.  The  heat  loss  is  taken  proportional  to  the  kinetic  energy  of  the  ejected  mass 
in  a 1:3  ratio. 


5.  The  ratio  of  specific  heats,  y,  is  taken  to  be  1.25. 

6.  The  multitube  catapult  can  be  treated  as  a constant  bore  catapult.  1'he  bore 
is  determined  by  that  tube  in  motion  at  the  moment  of  maximum  thrust. 


7.  The  retarding  rail  friction  force,  Ff,  is  proportional  to  the  square  of  the 
stroke,  i.  e.,  Ff  = K(X  - X0)2.  <2> 


DEFINITION  OF  SYMBOLS 


A = equivalent  bore  area 
B = burning  rate  coefficient^3^ 

C = weight  of  charge 

Cy  = mean  specific  heat  of  gas  at  constant  volume 
F = gas  thrust  (=  PA) 

Ff  = rail  friction  retarding  force 
F*  = impetus  of  propellant  (=  RTa)*3^ 


in. /spc/psi 
lb 

ft  - lb/lb/°F 
lb 

lk 

ft -lb  lb 


^ ^See  Appendix  I for  conversion  method. 

*^2^Vhile  this  is  8 typical  rail  friction  expression,  the  rail  friction  may  be  represented  by 
any  applicable  algebraic  function  of  X or  V. 

^^Theie  quantities  do  not  appear  directly,  but  their  product,  BF*  , appears  in  an  adjusted 
constant,  Q. 

^^These  are  nominal  quantities  that  are  used  in  the  derivations  of  F*  and  y , but  do  hot  appear 
.in  the  final  equations. 


g 


= gravitational  constant  (=  32.17) 

( 2 ) 

K = rail  friction  force  coefficient  where  Ff  = K(X  - X0)2 
L = equivalent  stroke 

M = ejected  mass 

N = weight  of  propellant  burned 

P = gas  pressure 

CBF* 

Q = a ballistic  parameter  (=  ) 

WA 

R = gas  constant  * 

T = temperature  of  gas  * 4 * 

Ta  = adiabatic  flame  temperature <4 5 

V = velocity  of  mass  with  respect  to  airplane,  along  path  of 

ejection 

V = acceleration  of  mass  with  respect  to  airplane,  along  path  of 

ejection 

VP/G  = acceleration  of  airplane  with  respect  to  earth 
VM/G  = acceleration  of  mass  with  respect  to  earth 
v = gas  volume  (=  12AX) 

v0  = initial  gas  volume 

W = web  thickness  for  a single  perforated  grain^5) 

X = equivalent  displacement  (=  travel  + X0) 

v0 

X0  = equivalent  initial  length  (= 

XB  = equivalent  muzzle  length  (=  X0  + L) 


ft/sec2 

lb/ft2 

ft 

s lugs 
lb 
ps  i 
fps 

ft- lb/lb°F 

°p 

^abs 
* ab  s 

fps 

ft/sec2 


ft/sec2 

ft/sec2 


in. 

ft 

ft 

ft 


Greek  Symbols 

a.  = load  factor  just  prior  to  ejection.  Also,  acceleration  in  g's  of  ejected  mass 
with  respect  to  airplane  due  to  inertial  and  gravitational  forces. 

( 2 \ 

While  this  is  a typical  rail  friction  expression,  the  rail  friction  may  be  represented  by 
any  applicable  algebraic  function  of  X or  V. 

(A  ) 

'These  are  nominal  quantities  that  are  used  in  the  derivations  of  F*  and  y,  but  do  not  appear 
in  the  final  equations. 

^S^The  equivalent  web  thickness  for  a 7*perforated  grain  is  taken  as  1.15  times  the  actual  web 
t hickne ss . 

^^By  convention,  a value  of  a greater  than  zero  represents  a net  force  aiding  the  gas  driving 
force,  F. 
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fi  = ratio  of  heat  transferred  to  catapult  to  kinetic  energy  of  ejected  mass 

y = ratio  of  specific  heats  ( = Cp/Cv ) 

y = a pseudo  ratio  of  specific  heats  such  that  y - 1 = (1  + (3)  (y  - 1) 

6 = angle  between  ejection  path  and  vertical  axis  at  time  of  ejection 

X.  v. 

= expansion  ratio  (=  — = — ) 

*o  vo 

Subscr i pts 

0 = initial  conditions  (defined  by  V0  a 0) 

p = peak  pressure  conditions 
b = "burnt"  conditions  (defined  by  N B C) 
m = "muzzle"  conditions 
bp  a black  powder 


METHOD  OF  SOLUTION 
Equations  for  the  System 

For  convenience  of  calculation,  the  airplane  is  taken  as  the  frame  of  reference 
and  all  forces  and  accelerations  are  summed  along  the  path  of  ejection. 

Figure  1 illustrates  the  forces  acting  upon  the  ejected  mass  when  the  ejection 
path  makes  an  angle  9 with  the  vertical  axis.  The  load  factor,  a,  is  the  ratio  of 
the  apparent  weight  to  the  actual  weight  of  the  ejected  mass  just  prior  to  ejection. 
This  apparent  weight  is  due  to  the  combined  effects  of  the  gravitational  and  inertial 
forces.  It  can  be  shown  that  a is  also  the  acceleration  in  g’s  of  the  ejected  mass 
with  respect  to  the  airplane  as  a result  of  these  inertial  and  gravitational  forces 
during  ejection. 

The  equation  of  motion  of  the  ejected  mass  with  respect  to  the  airplane  is  then 
F = MV  - aMg  + Ff  (1) 

A second  equation  of  motion  is 

V = X (2) 
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The  equation  of  motion  of  the  ejected  mass  with  respect  to  the  earth  is 


PA  - Ff  i Mg  COS  e = mv„  J7> 


Since  it  is  generally  true  that 


VM/G  = V +-VP/G 


We  may  find  VP/G  as  the  difference  between  equations  (1)  and  (3).  This  gives 


Vp/c  = Vm/g  ~ V = g(±  cos  d - a)(7) 


Therefore,  if  the  acceleration  of  the  mass  with  respect  to  the  airplane,  V,  is 
found  with  the  aid  of  the  computer,  the  acceleration  of  the  mass  with  respect  to  the 
earth,  VM/G,  can  be  found  from: 

VM/g  = V + g(±  cos  6 - a)<7>  /g) 


This  quantity  is  an  important  aeromedical  parameter. 


The  ballistic  equations  are: 
Equation  of  state 


Pv 

= FX  = NRT 

12 


Energy  equation 


NCVT.  = NCVT  + 1/2  MV2  + -^-MV2  - aMg(X  - X0)  + FfdX 


Equation  of  burning 


d fti  \ BP 
dtl  C / W 


Plu«  for  downward  ejection,  ainut  for  upward  ejection. 
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Neg.  #25758-1 
R- 1224 


VERTICAL  AXIS 

Figure  I.  Summation  of  forces  on  ejected  mass 


r 

t 
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Differentiation  of  equation  (7)  and  combination  with  equation  (2),  after  some  re- 
arrangement ,*8)  results  in: 


d N\  • - 

CF*  — - = XF  + yFV  + O'  - y)V(aMg  - Ff  ) 

dt  VC; 

This  may  be  combined  with  equation  (9)  to  give: 
QF  = XF  + yFV  + <y  - 7)V(iMg  - Ff) 

where 


Q = 


CBF* 

WA 


(10) 


(H) 


(12) 


It  may  be  noted  that  equation  (11)  is  valid  where  Ff  is  any  function  of  X. 


(8) 


See  Appendix  II  for  derivation  of  equation  10. 


* 

t 


The  solution  of  a specific  catapult  system  under  specific  operating  conditions 
can  be  achieved  if  six  constants  and  three  initial  conditions  are  known.  The  equations 
used  for  the  computer  are: 


F = MV  - aMg  + Ff 

(1) 

V = X 

(2) 

QF  = XF  + yFV  + (7  - 7)V(aMg  - Ff) 

(11) 

The  constants  for  the  system  are  M,  a,  K,  Q,  y,  and  y.  The  initial  conditions  for 
the  system  are  F0,  X0,  and  V0  = 0. 

Standard  difference  methods  are  used  with  a high-speed  digital  computer.  A 
Remington  Rand  model  No.  409-2A  computer  has  been  employed  for  this  problem  at  Frank - 
ford  Arsenal. 

The  output  of  the  conputer  tabulates  values  of  X,  V,  and  F for  each  time  incre- 
ment. (An  increment  of  five  milliseconds  has  been  found  to  be  satisfactory.)  The 
solution  is  halted  when  X attains  the  value  of  XB.  The  time  required  for  the  stroke 

and  the  velocity  at  muzzle,  V„,  may  be  noted.  The  peak  pressure,  Pp,  is  obtained  from 
the  greatest  value  of  thrust,  since  Pp  = Fp/A.  The  maximum  acceleration  of  the  mass, 

V,  may  be  obtained  from  the  greatest  difference  in  successive  values  of  velocity  since 
• • 

V ‘a  Av/At.  The  maximum  absolute  acceleration  of  the  mass,  VM/C,  may  then  be  found 
with  equation  (6). 


CONSTANTS  AND  INITIAL  CONDITIONS 

The  constants  and  initial  conditions  for  a given  case  may  be  evaluated  as  follows: 

1.  M is  the  ejected  mass  in  slugs. 

2.  a is  the  load  factor  just  prior  to  ejection  or  the  ratio  of  the  apparent 
weight  to  the  actual  weight,  Mg,  of  the  ejected  mass  just  prior  to  ejection.  Also,  a 
is  the  acceleration  in  g’s  of  the  ejected  mass  relative  to  the  airplane  as  a result  of 
the  inertial  and  gravitational  forces  during  ejection.  By  convention,  a positive  value 
of  a represents  an  acceleration  of  the  ejected  mass  relative  to  the  airplane  that  acts 
to  assist  the  ejection,  while  a negative  value  of  a.  represents  an  acceleration  of  the 
mass  relative  to  the  airplane  that  acts  to  retard  the  ejection. 

^^The  aethode  for  calculating  Xq  and  Xa  are  given  in  Appendix  I. 
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3.  K is  the  frictional  drag  coefficient  defined  by  the  relation: 


Ff  = K(X  - X0)2 


(13) 


The  value  of  K for  a given  system  flying  at  a given  altitude  at  a given 
speed  is  a function  of  the  aerodynamic  drag  on  the  ejected  mass,  the  guide  rail  design 
geometry,  and  the  coefficient  of  rail  friction.  For  a given  guide  rail  design  and 
coefficient  of  friction,  the  aerodynamic  drag  as  a function  of  stroke  may  be  determined. 
This  can  be  done  by  wind  tunnel  tests  or  aerodynamic  calculations.  The  frictional  drag, 
Ff,  may  then  be  determined  as  a function  of  stroke.  The  constant,  K,  is  then  chosen  by 
the  usual  means  to  give  the  best  fit  to  the  approx imate ly  parabolic  curve  of  Ff  vs  X. (10) 


4.  Q is  the  ballistic  constant  defined  by: 


0 = 


CBF* 

WA 


(12) 


The  values  of  C,  W,  and  A are  known  for  a given  system.  The  BF*  product 
is  a function  of  propellant  composition.  It  is  determined  from  a semi -empirical  cali- 
bration procedure  described  in  the  next  section. 


5.  y is  the  pseudo-ratio  of  specific  heats  defined  by  the  relationship 

y - i = (i  + /3)<y  - i)  (i4) 

y is  taken  as  1.25  for  propellant  combustion  products.  The  ratio  of  heat  loss  to 
kinetic  energy,  /3,  is  taken  as  1:3.  This  is  a common  assumption  for  guns  and  will  be 
used  until  experimental  data  become  available. 


6.  X0  is  the  initial  length  of  the  equivalent  catapult.  It  is  defined  as 


*0 


v0 

12A 


(15) 


7.  F0  is  the  initial  gas  force  at  the  instant  of  motion.  It  is  found  in 
either  of  two  ways: 


(a)  For  a < 0,  the  initial  force  may  be  taken  as  -aMg. 


(10)There  nay  be  cases  where  the  friction  function  is  best  fitted  by  an  expression  other  than 
the  commonly  used  parabolic  one.  The  expression  that  best  fits  the  data  should  be  used  (see 
Footnote  (2 ) ) . 


(b)  For  a ^ 0,  the  initial  force  may  be  taken  as  the  force  developed  by 
burning  the  booster  charge  of  black  powder  instantaneously  in  the  initial  volume. 
Then: ^ 1 1 ^ 


F„  = PoA 


12CbpF»  = bp 

v0  X0 


(16) 


where  F*bp  is  taken  as  80,000  ft- lb/lb. 


CALIBRATION  PROCEDURE 


The  purpose  of  the  calibration  procedure  is  to  determine  the  BF*  product  for  a 
particular  composition  of  propellant  when  burned  in  a propellent  actuated  device. 

A semi -empirical  technique  is  employed.  Firing  data  for  a standard  catapult  are  used. 
The  solution  for  that  catapult  is  generated  on  the  computer  for  an  assumed  value  of 
Q,  and  the  muzzle  velocity  thus  obtained,  V„,  is  noted.  Several  tries  of  Q are  made 
until  Vm  agrees  with  the  experimental  data.  The  value  of  Q for  this  trial  is  taken  as 
the  proper  value  of  Q.  As  a further  check,  values  of  maximum  thrust  and  time  of  stroke 
should  approximately  agree  with  experimental  values. {13) 

The  experimental  quantities  C,  W,  and  A are  known,  and  Q has  been  obtained  from 
a computer  match  of  muzzle  velocity.  The  BF*  product  may  therefore  be  found  from  the 
re lationship: 


CHF* 


The  value  of  BF*  is  a property  of  the  propellant  composition  and  may  be  used  to 
help  evaluate  Q for  a new  catapult  that  will  use  this  "calibrated"  propellant. 


If  the  BF*  product  is  determined  in  the  above  manner  from  two  or  more  catapults 
that  use  the  same  composition  propellant,  then  a mean  value  of  the  BF*  products  should 
be  taken  as  the  most  reliable  value. 

*"*Futur«  experimental  work  may  ahow  that  the  initial  gaa  force,  F0,  for  jero  or  poaitive  "g" 
firinga  ia  more  cloaely  dependent  upon  the  cartridge  caae  rupture  preaaure  than  the  black  powder 
charge  weight. 

In  any  caae,  it  nay  be  adviaable  to  ahlft  the  point  t = 0 for  caaea  of  large  poaitive  a from 
the  time  where  V = 0 to  the  time  where  the  cartridge  caae  rupturea.  Since  there  ia  a delay  time, 
tj,  from  when  the  latchea  are  releaaed  to  when  the  cartridge  caae  rupturea,  the  original  volume 
will  have  increaaed  and  the  tube  velocity  will  be  greater  than  zero  due  to  the  accelerating  field 
ag.  For  t = 0 at  cartridge  caae  rupture,  we  have  V0  = agtj  end  new  v0  = old  v0  + where 
Av  = )4  adgtd2.  Future  experimental  work  la  neceaaary  to  aacertain  t<j . 

M 9 1 

' 'This  will  probably  differ  from  the  BF * product  of  the  propellant  when  burned  in  a closed 
boab  or  vented  vessel. 

( i3  \ 

v ‘'Computer  calibrations  and  computer  curves  for  several  standardized  catapults,  together  with 
all  important  constants,  are  presented  in  Appendix  III. 
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EXTENSION  OF  SOLUTION 


After  the  performance  of  a proposed  catapult  has  been  evaluated  with  the  conqjuver 
and  found  to  be  satisfactory  for  the  normal  operating  conditions  (i.  e.,  a = -1,  K = 0), 
it  is  very  helpful  to  study  the  catapult  performance  under  combinations  of  g forces  and 
frictional  forces  that  may  be  expected  in  service.  Changes  in  the  ejection  velocity, 

V„,  the  maximum  acceleration,  V max;  the  peak  pressure,  Pp;  and  the  ejection  time,  tB, 
are  observed.  This  knowledge  is  important  for  rapid,  inexpensive  development  of  a new 
catapult . 


The  solutions  obtained  from  the  computer  analysis  all  assume  that  the  web  size  is 
sufficient  to  insure  burning  until  tube  separation.  This  is  so  because  the  term 


*s  eliminated  between  two  ballistic  equations,  (5)  and  (6),  and  is  never  entered 
into  the  computer  equations.  Hence,  the  computer  has  no  way  of  checking  when  the 
charge  has  been  consumed  (N  = C). 


If  the  web  thickness  is  sufficiently  large,  this  does  not  pose  a problem.  On  the 
other  hand,  if  the  web  is  of  just  sufficient  size  to  support  burning  until  tube  sepa- 
ration for  a normal  ejection,  the  grains  may  be  consumed  at  some  time  prior  to  tube 
separation  in  cases  of  large  g forces  and  large  frictional  forces  retarding  the  motion 
This  is  because  the  burning  rate  is  proportional  to  the  pressure,  and  the  pressure  (or 
force,  F)  is  higher  when  large  external  forces  retard  the  motion.  When  this  condition 
occurs,  an  estimation  of  muzzle  velocity,  VB,  with  "burnout"  may  be  made. 


The  burning  law  may  be  written  in  the  form 


(17) 


This  can  be  integrated  from  t = 0 to  t = tb  (where  N/C  =1)  and  rearranged  to  give 


t = tt 


/ 


CF* 

Fdt  = (1  - N0/C)  — 


t=0 


(18) 


Hence,  the  condition  of  "burnt"  can  be  expressed  by  equating  an  area  under  the 
F-t  curve  to  a constant. 
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The  term  N0/C  in  equation  (18)  is  a measure  of  the  charge  burned  at  the  instant 
of  motion.  This  will  increase  for  larger  negative  (retarding)  values  of  a.  and  can  be 
found  from 


N0  _aMgX0  CbpF*bP 


CF* 


CF* 


(19) 


The  second  term  on  the  right  side  of  equation  (19)  accounts  for  the  black  powder’s 
contribution  to  F0. 

In  order  to  use  equation  (18),  one  must  make  some  reasonable  assumption  from 
available  data  about  either  B,  F*.  or  the  point  of  "burnt."  For  the  solutions  shown 
in  Appendix  III,  it  was  assumed  that  either  VB  = Vb  when  a = -1  or  a reasonable  value 
of  F*  was  selected.  The  resulting  calculated  points  of  "burnt"  are  shown  on  the 
curves  in  Appendix  III  where  broken  lines  start  at  "burnt"  and  continue  to  the  "muzzle" 
point.  These  broken  lines  are  the  solutions  that  would  have  been  obtained  if  the  burn- 
out did  not  occur  until  the  "muzzle"  point.  A stepwise  numerical  solution  of  the  M1A1 
catapult  from  "burnt"  to  "muzzle"  is  shown  for  the  case  of  a = -7  to  illustrate  the 
expected  real  performance  when  V„>Vb. 

The  solution  from  "burnt"  to  "muzzle"  is  described  by  a pseudo-adiabatic  expansion. 

Pbvby  - P .v.7  (20a) 


FfaX^  = FaXj 


(20b) 


Fb  is  found  from  this  relationship  since  the  other  quantities  are  known.  The  work 
done  by  the  gas  during  this  expansion  is: 


F nX,  - FbXb 

i - y 


(21) 


This  energy  is  used  to  overcome  the  rail  friction  and  to  accelerate  the  catapult 
according  to  the  following  relationship: 


E = J F fdX  - 0Mg(X„  - Xb ) + 1/2  M(V.2  - Vb2) 
X-X„ 


(22) 


For  a known  friction  function,  this  equation  can  be  solved  for  Vm. 
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CONVERSION  FROM  ACTUAL  TO  EQUIVALENT  CATAPULT 


The  actual  catapult  of  several  concentric  tubes  is  reduced  to  an  equivalent  cata- 
pult of  constant  bore  for  convenient  analysis.  The  reference  bore  is  determined  by 
that  tube  in  motion .at  the  moment  of  maximum  thrust. 


Neg.  #25758-2 
R-1221 


ACTUAL  CATAPULT 
SYSTEM 


SYSTEM 


Figure  2.  Diagrammatic  representation  of  actual  and  equivalent 

catapult  systems 
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Figure  2 schematically  represents  an  actual  catapult  system  and  its  ballistic 
equivalent  for  the  case  of  maximum  thrust  occurring  during  motion  of  tube  1.  The 
following  symbols  are  used: 

V0  = initial  volume  of  actual  catapult  (including  the  annulus  volumes)  in 

A ! = effective  bore  area  of  tube  1.  This  is  determined  by  the  outside  in 

diameter  of  tube  1 at  its  center. 

Aj  = effective  bore  area  of  tube  2,  determined  in  a manner  similar  to  Aj.  in 

A = bore  area  of  equivalent  catapult  (shcnm  here  as  At).  in 


= stroke  length  of  tube  1 
j£2  = stroke  length  of  tube  2 
L = stroke  length  of  equivalent  catapult 
X0  = initial  length  of  equivalent  catapult 
XB  = total  length  of  equivalent  catapult  ( = X0  + L) 


For  A = A, : 

Aj 

L = + - — £ 2 

A1 


X.  = X0  + L 


For  A = Aj! 


L - — + -£2 

A 1 


X.  = X0  + L 


These  equations  may  be  extended  for  three  or  more  tubes. 
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DERIVATION  OF  EQUATION  (10) 


Pv 

12 

NCVT. 


FX  = I«T 

X 

NCVT  + 1/2  MV2  + - MV2  - aMg(X  - X0)  + / FfdX 
2 X = X0 


Differentiating  equation  (7)  with  respect  to  time,  we  have 
XF  + FV  = CRT. 


[n  _d_  [T\  T_  d_  ft$ 

r‘  Lcdt  [tT/  +t7^1cJ 


Rearrangement  of  equation  (8)  results  in 


Since 


_T  _ ( y - 1)  MV2  _ (y  - 1) /FfdX  (7  - l)aMg  (X  - X0) 

T.  ~ 2NRT,  NRT.  PRT. 


y - 1 


And 


y - i - (i  + 0X7  - i) 


Substitution  of  T/T.  from  equation  (B)  into  equation  (A)  gives 

<y  ~ 1)  MV2  (y  - 1)  f FrdX 


XF  + FV  = CF* 


/'N  _d_ 

lc  dt 


1 - 


2NF* 


NF* 


+ (y  - i)QMg(X  - X0)  II  (y-nr  (y 

NF*  I I 2NF* 


(y  - i)oMg(x  - x0> 

NF* 


Where 


RT.  = F* 


(7) 

(8) 

(A) 

(B) 

(C) 

(D) 


1 ) I FfdX 
NF* 

(E) 

(F) 
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When  the  indicated  differentiation  is  performed,  equation  (E)  takes  the  forn 


J FfdX  - N jFfdX 


Where 


1)  aMg  / 

F*  V 


NV  - (X  - X0)  N 


+ f 1 _ (T  - j)  . (7  - 1)  jf  FfdX  + (y  - 1)  aMg  (X  - X0)  1 
\ 2NF*  NF*  NF*  J 


d r d r dX 


Expansion  and  reduction  of  (G)  results  in 


XF  + FV  = - (y  - 1)  MW  - (y  - 1)  F,V  + (y  - 1)  aAfeV  + NF* 


Substituting  for  V in  equation  (I)  according  to  equation  (1),  i.  e.; 


V = — (F  + aMg  - Ff),  simplifies  equation  (I)  to 
M 


NF**XF+>FV  + (y-l)V  (aMg  - Ff ) - (y  - 1 ) V (aMg  - Ff ) 


Which  can  be  further  simplified  to 


CF*  — J = XF+yFV  + (y-y)V  (aMg  - Ff ) 


This  is  the  required  relationship. 


Computer  curves  for  several  catapults  with  g loads  that  retard  the  motion  are 
shown  in  the  attached  graphs  (Figures  3 to  9,  incl).  The  constants  used  for  the  calcu- 
lations are  summarized  in  the  following  chart. 


Propellant 
composit ion 


Unit 

MIAI  and 

M5  M2AI 

M3 

(M37  cartr idge ) 

TIP 

m 

lb 

312.0 

311.8 

365.2 

324.0 

363.4 

300.0 

-1.  -3, 

-1.  -3, 

-1.  -3, 

-1,  -3. 

-1.  -3. 

-1.  -3 

-5.  -7 

-5.  -7 

-5.  -7 

-5.  -7 

-5.  -7 

•5.  -7 

lb/ft2 

0 

0 

0 

0 

0 

0 

fps 

54.0 

43.2 

110.4 

44.0 

80.0 

26.8 

1.33 

1.33 

1.33 

1.33 

1.33 

1.33 

lb 

aMg 

aMg 

aMg 

aMg 

aMg 

aMg 

ft 

2.273 

2.160 

7.763 

3.296 

3.929 

1.130 

ft 

7.325 

6.493 

15.867 

7.531 

9.417 

2.883 

3.22 

3.005 

2.044 

2.285 

2.400 

2.55 

M6 

M6 

H8 

H8 

T8 

M2 

in/sec /ps i 

x ft- lb/lb 

105.3 

100.8 

66.8 

75.7 

110.0 

130.5 

lb 

0.1863 

0.1700 

0.441 

0.1214 

0.353 

0.0618 

in. 

0.1270 

0.1270 

0.110 

0.118 

0.130 

0.170 

in.  2 

2.86 

3.125 

2.426 

1.77 

3.733 

1.77 

Neg.  *25758-3 


Figure  3.  A comparison  of  computer  curves  with  experimental  curves  for  Catapult,  Aircraft  Personnel,  Ml  A I 


v in/ sec) 


X (FT)  VIFT/SEC)  F ( KIPS) 


Figure  6.  Computer  curves  of  thrust,  velocity,  and  displacement 
vs  time  for  Catapult,  Aircraft  Personnel,  M3 
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